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INTRODUCTION 
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i . 

This  report  has  two  objectives.  The  first  is  to  present 
some  new  analytic  solutions  for  source  region  EMP  from  a  sur¬ 
face  nuclear  burst  over  a  perfect  ground  in  the  quasi-static 
time  regime,  the  late-time  period  during  which  time  deriv¬ 
atives  in  Maxwell's  equations  can  be  neglected.  The  second 
objective  is  to  compare  the  existing  theory  of  source  region 
EMP,  as  embodied  in  these  analytic  solutions,  with  the  gross 
features  of  lightning  strokes  observed  on  the  Mike  shot 
(Ref.  1) . 

Quasi-static  source  region  EMP  has  been  investigated 
previously  by  Longmire  (Ref.  2),  Hill  (Ref.  3)  and  Wyatt 
(Ref.  4),  all  of  whom  noted  that  the  radial  component  of  the 
electric  field  would  be  smaller  than  the  polar  component,  and 
therefore  assumed  that  it  could  be  neglected  in  the  calcula¬ 
tions.  In  the  present  work  no  such  approximation  is  made.  The 
electrostatic  potential  is  expanded  in  spherical  harmonics 
(Legendre  polynominals)  to  obtain  a  set  of  ordinary  differ¬ 
ential  equations  for  the  radial  functions  appearing  in  the 
expansion.  Simple  air  conductivity  models  are  introduced  for 
the  early-time  regime  when  the  conductivity  is  dominated  by 
electrons  (Model  1) ,  the  late-time  regime  when  ions  dominate 
(Model  2) ,  and  for  the  intermediate  time  regime  when  electrons 
dominate  the  conductivity  close  to  the  burst  but  ions  dominate 
farther  away  (Model  3).  Using  these  simple  models,  analytic 
solutions  are  obtained  for  the  required  radial  functions  .  The 
results  are  then  evaluated  numerically  using  a  programmable 
calculator. 

From  the  present  work  it  is  found  that  neglect  of  the 
radial  component  of  the  electric  field  is  not  justified  when 
electrons  contribute  significantly  to  the  air  conductivity. 

In  addition,  it  is  found  (from  Model  2)  that  even  in  the  special 
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case  where  ions  dominate  the  air  conductivity,  the  previous 
approximate  solutions  obtained  by  the  above-referenced  authors 
are  only  qualitatively  correct.  The  results  from  Model  1  and 
Model  3  provide  new  insight  into  the  electromagnetic  field 
structure  when  electron  conductivity  cannot  be  neglected. 

Concerning  the  application  to  nuclear  lightning,  it  has 
previously  been  asserted  (Refs.  5—8)  that  the  gross  character¬ 
istics  of  lightning  observed  on  the  Mike  shot  are  clearly 
inconsistent  with  the  existing  theory  of  source  region  EMP . 
This  assertion  has  been  based  on  the  assumptions  that  (1)  the 
lightning  requires  fields  greater  than  100  kV/m,  whereas  the 
EMP  theory  predicts  fields  of  the  order  of  ~10-50  kV/m,  and 
that  (2)  the  observed  curvature  of  the  lightning  strokes 
cannot  be  explained  unless  ionic  air  conductivity  dominates 
over  electronic  conductivity  to  a  degree  significantly  greater 
than  the  existing  EMP  theory  would  allow.  Based  on  these 
assumptions,  it  has  been  concluded  by  others  that  the  observed 
lightning  provides  evidence  for  an  anomalously  rapid  electron 
attachment  mechanism,  which  lowers  electron  conductivity  and 
leads  to  higher  EMP  field  strengths. 

In  the  present  report  it  is  shown  that  these  assumptions 
do  not  appear  to  be  well  founded,  for  a  number  of  reasons. 
First,  it  is  pointed  out  that  existing  data  do  not  confirm 
the  hypothesized  100-kV/m  criterion  for  lightning  stroke 
initiation,  even  under  nonnuclear  conditions,  and  that  the 
presence  of  ionized,  conducting  air  may  also  facilitate  leader 
development  in  the  EMP  source  region.  Furthermore,  it  is 
found  that  the  observed  lightning  stroke  curvature  cannot  be 
shown  to  be  inconsistent  with  the  existing  EMP  theory  because 
of  (1)  the  hitherto  unrecognized  dominance  of  the  polar  compo¬ 
nent  of  the  electric  field  even  earlier  than  the  onset  of  ion 
conductivity  dominance,  (2)  the  continual  decline  in  the 
importance  of  electron  conductivity  during  the  period  of 
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lightning  stroke  evolution,  and  (3)  the  uncertainties  in  the 
present  theory  regarding  the  relative  contributions  of  elec¬ 
trons  and  ions  to  the  source  region  air  conductivity  at  late 
times.  Thus,  it  is  concluded  that  the  gross  features  of  the 
observed  lightning  may  very  well  be  consistent  with  the 
existing  theory  of  source  region  EMP .  Although  the  possibil¬ 
ity  of  some  anomalous  electron  attachment  mechanism  cannot  be 
ruled  out,  the  present  comparison  of  theory  and  observation 
does  not  find  evidence  in  support  of  its  existence. 
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II.  THE  MODELS 

For  times  later  than  about  10  ps  after  a  low-altitude 
nuclear  burst,  the  rate  of  change  of  the  Compton  current 
becomes  so  low  that  it  is  commonly  assumed  (Refs.  2-5)  that 
the  electromagnetic  fields  attain  a  quasi-static  equilibrium 
in  which  all  time  derivatives  in  Maxwell's  equations  can  be 
neglected  and  time  treated  as  a  parameter  rather  than  an 
independent  variable.  Then,  since 

"v“  x  ~E  =  -~B  ~  0  , 

the  electric  field  is  derivable  from  a  potential 

=  -"v<5>  • 

And  since 

\T*  (<jE  +  IT  )  =  -p  2:  0 

C 

(where  a  is  the  air  conductivity,  the  Compton  current  den¬ 
sity  and  p  the  charge  density),  the  potential  is  determined  by 
the  equation 

T-  (oV4>)  =  T- j“c  (i) 

subject  to  the  boundary  conditions  that  eT  be  finite  at  the 
origin  and  zero  at  infinity  and  (for  a  perfectly  conducting 
earth)  that  the  tangential  component  of  ~E  vanish  at  the  sur¬ 
face  of  the  earth. 
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In  the  present  work  we  will  assume  (Ref.  5)  that  the 
Compton  current  from  a  nuclear  surface  burst  may  be  approx¬ 
imately  represented  as 


Tc  =  rJoe~rA/r2  (2) 

where  r  is  a  unit  vector  in  the  radial  direction,  r  is  the 
distance  from  the  origin  and  X  an  effective  y-ray  mean  free 
path.  Monte  Carlo  studies  (Refs.  9,10)  indicate  that  the 
theta  component  of  will  generally  be  substantially  less 
than  the  radial  component,  except  at  r~X. 

The  air  conductivity  o  will  involve  contributions  from 
both  electrons  and  ions.  In  a  quasi-steady-state  situation 
the  electron  contribution  will  be  of  the  form  (Refs.  5,11) 

ae  ~  eueS/a  (3) 

where  e  is  the  electronic  charge,  u0  the  electron  mobility, 

S  the  local  ionization  rate  and  a  the  electron  attachment 
rate.  The  quasi-steady-state  form  for  ionic  conductivity 
will  be  (Ref.  5) 


cjj.  =  2e  Uj  (S/k)  (4) 

where  is  the  average  ion  mobility  and  k  the  ion-ion  recom¬ 

bination  rate  constant.  The  spatial  variation  of  o  =  og  +  Oj 
will  thus  be  determined  by  the  spatial  variation  of  the  ion¬ 
ization  rate,  which  we  will  assume  (Ref.  5)  to  be  of  the  form 

S  =  SQ  e-r//X/r2  .  (5) 
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In  addition,  the  electron  mobility  and  attachment  rate  will 
vary  weakly  with  the  local  value  of  the  electric  field 
strength  (Ref.  11).  However,  we  will  assume  that  this  varia¬ 
tion  can  be  neglected  in  solving  Eq .  (1),  provided  only  that 
the  values  of  u  and  a  used  in  the  final  analysis  are  consis- 
tent  with  the  calculated  average  local  electric-field  strength 
in  the  region  considered.  This  is  reasonable  since  the  ratio 
Ue/a,  in  Eq.  (3)  is  nearly  independent  of  electric-field  in  most 
situations  of  interest. 

Equations  (3-5)  imply  a  rather  complex  dependence  of  the 

-r/ A  2 

total  conductivity  on  range,  with  c  ~  e  /r  dominating  at 

e  -r /2  A 

short  ranges  (high  dose  rate)  and  Oj  ~  e  /r  dominating  at 

larger  ranges  (low  dose  rate) .  However,  at  sufficiently  early 
times  and/or  for  larger  yield  explosions,  the  source  strength 
(SQ)  will  be  large  enough  that  the  domination  of  oe  over  Oj 
will  continue  out  to  ranges  of  many  gamma  mean  free  paths. 

In  order  to  investigate  the  close-in  source  region  EMP  fields 
under  these  conditions  we  will  adopt  a  model  in  which  is 
ignored  altogether. 

Model  1  (Ml):  o(r)  =  aQ  e_r/^/r2  ,  (6) 

and  oq  will  scale  in  proportion  to  the  source  strength  SQ 
[defined  by  Eq.  (5) ] . 

Similarly,  at  sufficiently  late  times  and/or  for  lower 
yield  explosions,  the  source  strength  will  be  so  low  that 
will  dominate  over  a&  everywhere  except  within  a  few  hundred 
meters  of  the  detonation  point.  For  this  situation  we  will 
adopt  a  model  in  which  og  is  ignored  altogether. 

Model  2  (M2):  o(r)  =  aQ  e“r/2X/rx,  (7) 

and  ctq  will  scale  as  the  square  root  of  SQ. 
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Finally,  for  the  intermediate  situation  where  oQ  dominates 

out  to  r  1  km  or  so,  with  dominating  at  the  larger  ranges, 

we  note  that  since  X  is  typically  about  300  m,  o(r)  varies 

2 

approximately  as  1/r  at  the  shorter  ranges  (r  55  1  km)  and 

-r /2  A 

approximately  as  e  '  at  the  longer  ranges  (r  S  1500  m) . 

For  this  situation  we  adopt  a  simplified  conductivity  model 
which  exhibits  this  dominant  behavior: 

Model  3  (M3):  o(r)  =  aQe~r‘/2X/r2 .  (8) 

Here  the  value  of  oQ  must  be  obtained  by  a  fit  to  the  estimated 
true  form  of  c(r)  =  oe(r)+aI(r).  Such  a  fit  is  shown  in 
Figure  1  for  the  case  of  a  10-MT  surface  burst  at  t  ==1  ms.  The 
parameters  used  in  constructing  these  curves  will  be  discussed 
later. 
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1.  Comparison  of  Model  3  to  expected 
conductivity  for  10-MT  burst  at 


III.  THE  SOLUTIONS 


Equation  (1)  can  be  rewritten  in  spherical  coordinates  as 
(we  ignore  a  6 (r)  on  the  right-hand  side  which  does  not  effect 
the  solution  for  r  0) 


2  9r  V 


r2o (r ) 


3_ 

3r 


4>  (r ,  0)  + 


o  (r) 


r  sin  6 


i  «<r,9))  = 


-  J 


(9) 


The  angle  0  is  measured  from  the  vertical,  so  that  9=0 
directly  above  the  burst  and  0  =  r/2  at  the  surface  of  the 
earth.  For  0  ^  0  <  tt/2  the  9 -independent  right-hand  side  of 
Eq.  (9)  may  be  written  as  (Ref.  12) 


oo 

(J-T(r>)=  -  i  (Mf<r>)  2A2i+l  P2i 


2i+l (cDS  6) 


1=0 


where  the  P2i+i^^  are  Legendre  polynomials  and  (l  =  2i+l) 

.  _  (  1  \U-l)/2/(2i+l)  (1-2)  !  !  \ 

■l-\2}  ^  2(fTTp  ) 

If  we  also  represent  $  (r ,  Q)  as 

oo 

Mr,  0)  =  ^J^i+l (r)  P2i+1  (cOS  0) 

i=0 


(10) 


(11) 


Then  the  boundary  condition  on  Er  at  6  =  90°  is  met  automati¬ 
cally  and 


dr 


2  t>lir)  +  r  dr  h(r)  +  a( 


1_  (do  (r)\ 
(r)  \  dr  / 


d4>£  (r) 
dr 


£ (£+1) 


*£(r)  = 


Jk  J  (r) 
X  o  (r) 


:i2) 


From  the  defining  Eqs.  (2,6, 7, 8)  we  may  then  introduce 


FTFT  3r  a(r)  = 


and 


r-J(r)/o(r)  = 


(V0o> 


-a  *  f) 

(Ml) 

(13a) 

-(h + *) 

(M2) 

(13b) 

-/J_  +  2) 

\  2A  r  / 

(M3) 

(13c) 

(Ml) 

(14a) 

X  -r/2X 

F  e 

(M2) 

(14b) 

e"r/2X 


o'  o 


(M3) 


(14c) 


into  (12)  to  obtain  the  appropriate  radial  equation  for  each 
model . 
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For  Model  1  the  homogeneous  equation  corresponding  to 
Eq.  (12)  has  two  independent  solutions  (Ref.  13)  of  the  form 


tl(r)  =  \/  r/2  X  er/2X 

1  ^ +  1/2  ^r/2  ‘  ■' 

(15a) 

£  2  (r )  =  s[~t/2  X  er/2X 

KJ+l/2<r/2>) 

(15b) 

where  an<^  K£+l/2  ^  are  i-ed  spherical  Bessel 

functions.  The  solution  to  Eq.  (12)  is  then  given  by 


*4(r)  =  4XA^(Jo/ao)/^72X  e 

r/2X 


r/2  X 


cx> 

Iz+i/2(r/2X)fdyKi+i/2{y)Jy  e~y 

r/2  X 


+  Kl+l/2(r/2X)J  dy  JZ+l/2ly)  /y 


-y 


(16) 


2 

It  can  be  shown  from  Eq.  (16)  that  <f>.  (r)  behaves  like  r 

2  * 

[or  r  £n(r)  for  l  =  1]  as  r  -*+  0  and  like  r  as  r  oo.  Note 

that  this  result  (and  the  basic  model  as  well)  is  somewhat 
unphysical  in  that  the  electric  field  is  nonzero  at  r  -*•  °j. 

This  is  at  least  partially  a  result  of  the  unphysical  assump¬ 
tion  that  Oj  <<c  CTe  (so  that  r*3/cr  =  a  constant)  even  as  r  -*■  °°. 
However,  we  will  only  apply  Model  1  at  short  ranges  from  the 

burst,  where  a  >>  a_. 

e  I 

For  Model  2,  the  homogeneous  variant  of  Eq.  (12)  has  two 
independent  solutions  of  the  form  (Ref.  14) 


<f>£l(r)  =  r?  M(p,2p+l,r/2X)  (17a) 

<J,fc2^r)  =  rP  U(P»  2p+l,r/2\)  (17b) 
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where  p  =  J £  (£+1)  and  M(p, 2p+l,r/2A)  and  U (p, 2p+l , r/2 A ) 
are  confluent  hypergeometric  functions.  The  solution 
to  Eq.  (12)  may  then  be  written  as  [T  (p)  is  the  gamma 
function] 


4>j,(r)  XA«-(aQjpr  (2p)  (r/2A)P 

dy  U (p , 2p+l ,y ) y^e 

A 

r/2  A 

+  U  (p  ,2p+l ,  r/2A  )  J  dy  M  ( p,  2p+l  ,y )  ype  2y 

o 

From  Eq.  (18)  it  can  be  shown  that  ^(r)  behaves  like  r  as 
r  0  and  like  a  constant  as  r  -*-oo.  In  this  case  the 
electric  field  vanishes  at  r  =  oo. 

Finally,  for  Model  3,  the  homogeneous  equation  corres¬ 
ponding  to  Eq.  (12)  has  independent  solutions  (Ref.  13) 

4>u(r)  =  /^7TTer/4>'Iji  +  1/2  (r/4A) 

4>*2(r)  =  /r/4  A  er/4XKp+1/2(r/4A) 


X  |m(p ,2p+l,  r/2 A) 


/ 


(18) 


(19a) 


(19b) 
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The  solution  to  Eq.  (12)  is  then  given  by 


>£ (r)  =  16A^  A 


(a 


r/4 A  e 


r/4  A 


00 

X  I£+l/2(r/4X)  f  dy  KH+l/2{y)  e* 

r/4  A 


r/4  A 

+  Vl/2(r/4M  J  dy  *4+1/2  (y)  e'3Y  * 

0 

2  2 

In  this  case  ^(r)  behaves  like  r  (or  r  £n  r)  as  r  -►  0  and 
like  a  constant  as  r -►oo,  thus  the  electric  field  vanishes 
at  r-*~oo. 

The  Eqs.  (10,11,16,18,20)  provide  a  set  of  exact  closed  - 
form  solutions  for  each  of  our  three  models.  The  electric 
fields  can  then  be  obtained  from  the  equations 


E0  S1r  /Lj  ^i+l (r)  (d(cos  8)  P2i+l(cos  0)) 
i=0 


00 


Er  a?  ^2i+l(r  /  P2i+1 (cos  e>  * 


In  Eq.  (22)  we  must  have  expressions  for  d4>p(r)/dr  for  each 
model  (l  =  2i+l) .  These  can  be  found  easily  by  differentiating 
Eqs.  (16,  18,  20)  and  using  recursion  relations  (Refs.  13,  14) 
for  the  derivatives  of  the  special  functions. 
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IV.  PHYSICAL  RESULTS 


In  order  to  generate  specific  results  for  the  electric 
fields  predicted  by  each  model,  we  must  choose  numerical 
values  for  the  parameters  X,  JQ  and  c  .  The  value  of  oq  will, 
of  course,  be  different  for  each  of  the  three  models. 

Here  we  will  take  X  =  320  m  for  the  effective  gamma  dose 
attenuation  length.  For  the  Compton  current,  a  review  of  the 
Monte  Carlo  calculations  by  Sargis  et  al.  (Ref.  9)  indicates 
that  at  t  ~  1  ms  and  500  m  £  r  £  2900  m 

~  -)  _2 4  3 

r*Jc(A/m  )  2;  -5.6  x  10  S  (ion-pairs/m  /s) 

where  S  is  the  local  ionization  rate.  On  the  other  hand,  the 
Monte  Carlo  calculations  by  the  Mathematical  Applications  Group, 
Inc.  (MAGI)  (Ref.  10)  are  more  consistent  with  a  proportionality 

r-J^.(A/ni  )  ^  -1.2  x  10_  S  (ion-pairs/m  /s)  . 

For  the  present,  we  will  use  an  intermediate  value,  choosing 

r*Jc(A/m2)  =  -8.2  x  10  S (ion-pairs/m3/s) .  (23) 

And  the  uncertainty  associated  with  this  value  is  evidently 
at  least  +40  percent. 

In  the  electron  conductivity  [see  Eq.  (3) ]  we  will  use 
(F.ef .  11) 


-  r  -1 

a  =  1 . 5  x  10  s 
y e  =  0.25  m2/V-s  . 
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These  values  are  appropriate  for  electric  field  strengths  of  a 
few  tens  of  kilovolts  per  meter  (kV/m)  and  air  with  ^2  percent 
moisture  content.  They  are  each  thought  to  be  uncertain  by 
+30  to  40  percent.  From  Eq.  (3)  we  then  obtain  (the  electric 
field  dependence  of  the  ratio  Ve/a  is  not  important  in  slightly 
moist  air) 

oe(mho/m)  =  2.7  x  10  2a  S ( ion-pairs/m3/ s) .  (24) 

Using  these  numerical  values  in  Model  1  (electron-domin¬ 
ated  conductivity)  we  obtain 

Ml:  J  /a  -  -30  kV/m 
o  o 

a  value  which  must  be  regarded  as  uncertain  by  at  least  +50  per¬ 
cent.  With  this  value  of  JQ/ao  and  A  =  320  m,  the  calculated 
electric  fields  for  Model  1  are  shown  in  Figure  2.  The  radial 
field  Er  is  calculated  at  6  =  0°  and  the  polar  field  E0  =  90° 

(at  the  ground) ,  since  they  are  each  expected  to  have  their 
largest  values  at  these  angles. 

In  the  ionic  conductivity  [see  Eq.  (4)]  we  will  use  the 
values  (Refs.  5,15) 


k  =  2  x  10~12  m3/s 
=  2.6  x  10~4  m2/V-s 

and  we  roughly  estimate  the  uncertainties  of  these  values  to 
be  about  +30  percent  (for  Uj)  to  +50  percent  (for  k) .  From 
Eq.  (4)  we  have 

°I  (mho/m)  =  5.9  x  10  17  S  (ion-pairs/m3/s) 1/'2  .  (25) 
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320  m,  -J  /rj  =  30  kV/m) 


For  Model  2  (ion-dominated  conductivity )  the  ratio  (J  /o  ) 

1/2  °  ° 
will  be  proportional  to  So  ,  which  is  essentially  the  square 

root  of  the  source  strength  [see  Eq.  (5)].  From  a  review  of 

the  Monte  Carlo  results  (Refs.  9,10)  we  find  a  value  of  S 

26  ^ 

at  t  ^  30  ms  for  a  10-MT  burst,  assuming  2  x  10  source  neu- 

2  8 

trons  per  MT  (Refs.  16,  17),  of  (30  ms)  -  7.2  x  10  ion- 
pairs/m/s.  For  this  particular  case  we  then  obtain  a  value  of 
Jq/Oq  for  Model  3  (again  with  at  least  +50  percent  uncertainty) 

M2:  <J0/a0)  =  -116  kV/m. 

The  electric  fields  E  (0=0°)  and  Eo(0  =  90°)  calculated 

r  t) 

from  Model  2  with  this  value  of  J  /o  and  X  =  320  m  are  shown 

o  o 

in  Figure  3.  For  other  burst  yields  and/or  times  after  deton¬ 
ation,  the  appropriate  value  of  S  will  be  different  and  the 

1/2  ° 

fields  will  scale  as  S  ^ 

o 

It  is  also  important  to  compare  the  calculated  values  of 
ag(r)  and  aT(r)  at  some  fairly  close-in  range,  say  r  =  500  m, 

for  this  choice  of  S  .  We  find  0  (500  m)  c:  o_(500  m) 

—  5  O  S  I 

—  1.5  x  10  mho/m.  Clearly,  Model  2  would  be  of  question¬ 
able  validity  for  much  higher  values  of  S  .  Note,  however, 
that  for  substantially  lower  yields  Model  2  would  become 
applicable  at  much  earlier  times,  perhaps  as  early  as  ~30  ys 
for  a  20-KT  burst  (however,  the  displacement  current  will  also 
start  to  become  important  for  such  low  yields  and  early  times) . 
The  presence  of  a  highly  conducting  expanding  fireball  would 
also  detract  from  the  applicability  of  the  model  for  very 
large  yields.  For  a  10-MT  burst  at  t  ~  30  ms,  the  fireball 
will  have  reached  about  500  m. 

As  for  the  limits  of  applicability  of  Model  1  we  may 
note  from  the  Monte  Carlo  work  (Refs.  9,10)  that  the 
source  strength  at  t  ~ 30  ys  for  a  10-MT  burst  would  be  about 
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Figure  3.  Numerical  results  from  Model 
and  -E„(90°) (assumes  A  =  320 


2.5  x  10"^  ion-pairs/m/s.  For  this  value  of  SQ  the  electronic 
and  ionic  contributions  become  comparable  at  r  =  1700  m.  Thus, 
even  for  a  10-MT  burst.  Model  1  must  have  a  fairly  limited 
range  of  applicability  (i.e.,  t  £  30  us,  r  ~  1500  m) .  For 
lower  yields  (lower  S  )  the  requiremenc  that  og  dominate 
over  0^  will  be  even  more  restrictive. 

Model  3  is  intended  to  at  least  partially  bridge  the  gap 
between  Models  1  and  2.  As  shown  in  Figure  1,  this  model 
provides  a  good  fit  to  the  total  conductivity  from  a  10-MT 
burst  at  t  ~  1  its.  The  source  strength  used  in  Figure  1  was 
So  =  1.1  x  10 J  ion-pairs/m/s.  At  this  source  strength  the 
electron  and  ion  conductivities  become  comparable  at  about 
1  km.  Fitting  the  model  form  of  the  conductivity  to  the 
expected  true  form  gave 

M3:  3  /o  =  -63  kV/m  . 

o  o 

The  electric  fields  E  (0°)  and  EQ  (90°)  calculated  from 

r  y 

Model  3  with  this  value  of  ^Q/°0  are  displayed  in  Figure  4. 

In  order  to  examine  the  limits  of  validity  of  Model  3, 

we  have  shown  in  Figure  5  a  fit  of  the  model  form  of  the 

conductivity  to  the  exoected  true  form  for  S  =  1.6  x  10^ 

29  ° 

ion/pairs/m/s  and  for  SQ  =  3.6  x  10  ion-pairs/m/s,  appropri¬ 
ate  for  a  10-MT  burst  at  t  ~  300  us  and  t  ~  3  ms,  respec¬ 
tively.  It  may  be  seen  that  the  model  still  provides  a  fair 
fit,  differing  from  the  desired  form  by  less  than  50  percent 
over  a  range  of  five  orders  of  magnitude  in  conductivity.  Thus, 
for  a  10-MT  burst,  we  might  expect  the  electric  fields  to 
exhibit  a  qualitative  similarity  to  those  shown  in  Figure  5 
for  the  interval  300  us  S.  t  S,  3  ms.  Also,  since  the  value  of 

a  used  for  these  fits  decreases  by  about  a  factor  of  3.3  from 
o 

t  —  300  us  to  t  a  3  ms,  while  SQ  (and  JQ)  decreases  by  a  factor 
of  4.4,  the  EMP  fields  (proportional  to  J  /aQ)  are  changing 
very  slowly  in  magnitude. 


<0  =  0°)  and  -Eg  (6  =  90°)  (kV/m) 


r  (m) 

Figure  4.  Results  from  Model  3  for  Er  and  -Eg  (assumes 
X  =  320  m,  -J  /a  =63  kV/m) . 


r 


Figure  5.  Comparison  of  Model  3  to  expected  true 
form  for  a  10-MT  burst  at  t  ~  300  us 
and  at  t  ~  3  ms. 
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Finally,  we  must  address  the  limits  of  applicability  of 

the  quasi-static  hypothesis,  upon  which  all  of  the  present 

calculations  have  been  based.  There  are  two  conditions.  The 

displacement  current  must  be  much  less  than  the  conduction 

current;  e.g.,  o/e  >>  t  ^ .  Also,  the  electric  field  must  be 

derivable  from  a  scalar  potential.  The  first  condition  is 

more  easily  met  for  large  yields  at  short  ranges  and  late  times. 

It  becomes  marginal  for  a  10-MT  burst  at  t  i  3  us  and  r  Z  3  km. 

For  a  20-KT  burst,  the  condition  becomes  marginal  at  t  £.30  us 

and  r  Z  2  km.  The  second  condition  can  be  shown  to  reduce  to 

2 

a  requirement  that  a  «  t/uQ^  ,  where  uD  is  the  permeability  of 
free  space.  This  condition  is  more  easily  met  for  lower  yields 
at  later  times  and  larger  ranges.  For  a  10-MT  burst  it  becomes 
marginal  at  t  £,  100  us  and  r  £  1  km.  For  a  20-KT  burst  it 
becomes  marginal  at  t  £.  100  us  and  r  £.  500  m.  It  may  be  seen 
that  this  last  condition,  combined  with  the  requirement  that 
ae  >>  Oj,  appears  to  place  a  very  severe  restriction  on  the 
applicability  of  Ml. 
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V.  COMPARISON  WITH  OTHER  WORK 

Longmire  (Ref.  2)  was  apparently  the  first  to  address  the 
problem  of  quasi-static  source  region  EMP.  He  developed  an 
approximate  solution  to  Model  2  by  noting  that  Er  should  be 
very  small,  also  changing  sign  at  r  ~  2X.  On  this  basis,  an 
explicit  expression  was  obtained  for  the  azimuthal  magnetic- 
field  at  0  =  90°. 


B0  =  r(*'Tc{t))  ■ 

In  analogous  developments,  Hill  (Ref.  3)  and  later  Wyatt 

(Ref.  4)  also  gave  approximate  expressions  for  the  polar 

electric  field.  Wyatt's  result,  including  a  polar  component 

of  J  ,  is 
c 


This  approximation  is  also  plotted  in  Figure  3  as  the  dashed 
line.  It  may  be  seen  that  the  approximate  result  is  fairly 
reasonable,  particularly  in  those  regions  where  Er  and  dEr/dr 
are  smallest. 

Hill  applied  this  result  to  estimate  E^  for  a  10-MT  burst, 
at  t  ~  3  ns,  but  assumed  a  substantially  lower  dose  rate  than 
would  be  indicated  by  the  Monte  Carlo  work  (e.g.,  Hill  uses 
SQ  ~  2.4  x  102*  ion-pairs/m/s,  compared  to  Sq  ~  3.6  x  1029  ion- 
pairs/m/s  from  the  Monte  Carlo  work.)  Because  of  this 
assumed  dose  rate,  the  ion  conductivity  was  dominant  in  Hill's 
problem,  as  required  for  Model  2  to  apply. 
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In  Wyatt's  application  of  the  approximate  solution  to 

29 

Model  2,  a  source  strength  of  SQ  =  7.2  x  10  ion-pairs/m/s 

was  assumed  for  a  10-MT  burst  at  t  ~  10  ms.  This  dose  rate  is 
about  3  to  4  times  larger  than  would  be  indicated  by  the  Monte 
Carlo  work.  At  this  dose  rate  Model  3,  with  both  ionic  and 
electronic  contributions  to  the  air  conductivity,  should  have 
been  used  rather  than  the  approximate  solution  to  Model  2. 
Also,  it  appears  that  Wyatt  miscalculated  the  ionic  conduc¬ 
tivity  by  a  factor  of  2,  apparently  not  including  the  con¬ 
tributions  of  both  positive  and  negative  ions,  and  used  an 
anomalously  large  ion-ion  recombination  rate  for  clustered 
ions. 


VI.  APPLICATION  TO  NUCLEAR  LIGHTNING 

This  work  was  motivated  primarily  by  the  recent  interest 
(Refs.  3-8)  in  lightning  strokes  triggered  by  nuclear  explo¬ 
sions  (Ref.  1).  On  the  Mike  shot,  a  10-MT  surface  burst, 
five  lightning  strokes  were  observed  at  ranges  from  900  m 
to  1380  m  from  ground  zero.  The  lightning  strokes  began  near 
the  ground,  apparently  at  t<l  ms,  and  travelled  upward  follow 
ing  roughly  concentric  paths  about  the  burst  point  throughout 
the  time  frame  of  1  ms£t  ^10  ms.  It  is  thought  that  they  may 
have  been  initiated  at  the  tips  of  communications  antennas  or 
other  structures  at  varying  distances  from  ground  zero.  A 
detailed  drawing  (Ref.  1)  showing  the  progression  of  the 
brightest  stroke  (at  r  =  900  m)  is  given  in  Figure  6.  In 
every  case,  the  cloud  level  was  low  enough  that  the  strokes 
could  not  be  observed  at  elevations  higher  than  30°  to  45° 
above  the  horizontal,  as  viewed  from  the  burst  point.  Also, 
a  photograph  taken  at  36  ms  indicates  (Ref.  1)  that  only 
one  of  the  strokes  actually  reached  cloud  level. 

From  a  study  of  several  examples  of  triggered  ground-to- 
cloud  lightning  under  nonnuclear  conditions.  Pierce  (Ref.  18) 
proposed  an  empirical  rule  that  EQh  =  (ambient  field)  x 
(structure  height)  >  1  MV.  Such  a  scaling  of  field  strength 
vs  structure  height  does  not  seem  implausible,  since  the  poten 
tial  drop  within  a  few  structure  radii  of  the  tip  .is  also  of 
the  order  of  EQh,  and  this  potential  drop  should  be  important 
in  initiating  any  discharge.  Pierce's  rule,  with  an  assumed 
structure  height  (h)  of  ^10  m,  led  Uman  et  al.  (Ref.  1)  to 
suggest  that  a  100-kV/m  EMP  field  would  be  needed  to  explain 
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the  Mike  lightning  (assuming  that  the  radiation-induced  ioniza¬ 
tion  is  unimportant).  However,  Pierce's  data  set  did  not  include 
any  relevant  examples  of  triggered  lightning  with  h  £  ICO  m. 

(One  example  at  h  =  70  m  apparently  had  no  real  data  on  Eq. 

Two  other  examples  of  lightning  strikes  to  in-flight  aircraft 
having  dimensions  of  5  to  36  m  gave  Eh—  36-75  kV,  a  factor 
of  13  to  28  lower  than  the  proposed  1-MV  criterion!)  The 
remaining  examples  generally  appear  to  conform  to  Pierce's 
suggested  rule,  but  show  a  scatter  of  about  a  factor  of  3  in 
EQh,  while  Eq  ranges  only  from  3  kV/m  to  20  kV/m  and  h  varies 
only  from  100  m  to  400  m. 

Overall,  it  cannot  be  said  that  the  data  demonstrate  the 
validity  of  a  scaling  law  of  the  form  EQh  —  constant  for  trig¬ 
gered  lightning,  or  that  they  can  be  extrapolated  with  better 
than  order-of-magnitude  accuracy  to  estimate  the  value  of  Eq 
needed  to  trigger  lightning  from  a  structure  with  h  r  10  m. 
Furthermore,  after  having  examined  a  variety  of  empirical 
schemes  for  plotting  and  extrapolating  the  data  cited  by  Pierce, 
this  author  holds  the  view  that  they  are  not  inconsistent  with 
a  required  field  strength  anywhere  in  the  general  range 
15  kV/m  &  Eq  &  150  kV/m  for  triggered  lightning  from  a  struc¬ 
ture  with  h  —  10  m.  The  EMP  fields  which  we  have  just  calcu¬ 
lated  are  certainly  within  this  range,  particularly  when  we 
recall  that  the  values  of  J0/°0  employed  are  probably  uncer¬ 
tain  by  at  least  +50  percent. 

It  clearly  would  be  very  useful  to  have  data  on  the  pres¬ 
ence  or  absence  of  upward  propagating  strokes  from  structures 
with  h  $  10  m  in  ambient  fields  of  the  order  of  20  kV/m  to 
100  kV/m.  However,  field  strengths  at  the  earth's  surface 
beneath  storm  clouds  do  not  frequently  exceed  20  kV/m  (Ref.  19) , 
and  triggered  ground-to-cloud  strokes  from  structures  with 


h  ~  10  m  are  apparently  also  very  uncommon.  Thus,  it  may  be 
difficult  to  find  data  in  the  regime  of  interest  from  observa¬ 
tions  of  naturally  occurring  ground-to-cloud  lightning  (the 
absence  of  massive  nuclear  radiation  would  also  be  a  concern) . 

Stronger  surface  fields  do  arise  when  a  cloud-to-ground 
stroke  approaches  the  earth,  and  such  fields  very  commonly  trig¬ 
ger  upward-propagating  connecting  leaders  from  short  (h  &  10  m) 
grounded  conductors.  This  is,  in  fact,  the  protective  mechanism 
of  the  lightning  rod.  It  is  generally  believed  (Refs.  19,20) 
that  such  upward  strokes  will  begin  when  the  enhanced  field  at 
the  conductor  tip  reaches  about  3  MV/m.  In  an  ambient  field 
of  30  kV/m,  this  would  require  a  prolate  spheroidal  structure 
(Ref.  21)  with  an  aspect  ratio  of  h/r  cr  20. 

Although  it  is  quite  plausible  that  the  enhanced  fields  at 
the  tips  of  projecting  conductors  within  1  to  2  km  from  a  10-MT 
burst  at  t  $  3  ms  will  exceed  3  MV/m,  so  that  some  type  of  dis¬ 
charge  phenomenon  should  be  expected,  it  is  not  clear  whether 
an  upward  propagating  streamer  will  result,  as  in  the  lightning 
rod  phenomenon.  The  applied  field  from  a  downward  propagating 
cloud- to-ground  stroke  is  rapidly  increasing  in  time,  while 
that  from  an  EMP  environment  is  nearly  constant.  It  is  also 
known  from  laboratory  experiments  (Ref.  22)  that  steady  fields 
only  slightly  above  breakdown  over  small  volumes  near  the  tip 
of  a  positive  electrode  generally  tend  to  produce  intermittent 
coronas  rather  than  leader  channels,  with  the  corona- to- leader 
transition  being  inhibited  by  space  charge  fields  which  oppose 
the  applied  field  in  the  region  between  the  space  charge  and 
the  electrode.  In  the  lightning  rod  phenomenon,  such  space 
charge  fields  are  probably  overcome  by  the  continual  increase 
in  applied  field  as  the  downward  propagating  stroke  approaches 
the  earth.  However,  in  the  nuclear  environment  the  space 
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charges  may  be  reduced  by  conduction  currents  through  the  sur¬ 
rounding  ionized  air,  a  possibility  to  which  we  will  shortly 
return. 

Before  addressing  the  potential  impact  of  ionization  and 
air  conductivity,  there  is  another  source  of  potentially  rele¬ 
vant  information  which  should  be  mentioned.  Data  obtained  from 
laboratory  studies  of  the  breakdown  of  long  (up  to  28  m)  posi- 
tive-point-to-plane  gaps  in  air  (Ref.  22)  have  led  to  the 
adoption  of  empirical  rules  of  the  form  VD (kV)  500  d^'5  for 
the  voltage  needed  to  break  down  a  gap  of  length  d(m).  Taking 
d  ~  600  m,  about  the  length  of  the  longest  arcs  in  the  Mike 
shot,  one  obtains  V0  —  12  MV,  whereas  the  present  EMP  calcula¬ 
tions  give  a  potential  difference  of  about  10  MV  (with  a 
+50  percent  uncertainty)  over  this  path. 

Finally,  we  should  consider  the  potential  effect  of  ioniza¬ 
tion  and  air  conductivity.  As  mentioned  previously,  a  leader 
channel  from  a  positive  point  electrode  is  typically  preceded 
by  an  "impulse  corona"  which  creates  a  positive  space  charge. 

The  effect  of  this  space  charge  is  to  suppress  further  elec¬ 
tron  avalanching  until  (Ref.  23)  the  energy  deposited  at  its 
stem  by  the  electron  avalanche  process  can  be  thermalized  to 
bring  the  gas  to  a  temperature  of  ~  3000°K;  at  this  point,  thermal 
detachment  of  electrons  previously  attached  to  oxygen  molecules 
raises  the  corona  stem  conductivity  to  the  point  where  the  space 
charges  can  begin  to  relax  and  also  allows  thermal  expansion 
of  the  gas,  which  increases  the  E/N  ratio  allowing  further  ava¬ 
lanche  development.  A  leader  channel  then  develops  from  the 
corona  stem  (provided  that  the  voltage  has  not  fallen) ,  and 
this  leader  channel  propagates  away  from  the  conducting  point, 
preceded  by  a  new  "leader  corona,"  which  feeds  current  into 


35 


the  advancing  leader  channel.  Depending  on  the  applied  volt¬ 
age,  the  voltage  drop  along  the  leader,  and  the  space  charge 
fields  of  the  corona,  the  leader  may  progress  to  bridge  the 
entire  gap,  or  it  may  eventually  stop  in  midair. 

According  to  this  picture  of  discharge  development,  the 
first  hurdle  to  overcome  in  creating  a  leader  channel  is  to 
relax  the  space  charge  fields  from  the  impulse  corona.  In 
the  nuclear  environment  there  is  a  new  mechanism  for  achieving 
this.  Since  the  air  conductivity  at  relevant  times  and  ranges 
is  of  the  order  of  3  x  10  mho/m,  the  space  charge  should  be 
relaxed  by  conduction  currents  in  times  of  the  order  of  a  few 
times  e^/o  —  0.3  ys.  Thus,  one  would  expect  leader  channel 

initiation  to  be  very  different  and  possibly  easier  in  the 

nuclear  environment  than  in  un-ionized  air. 

Similar  considerations  should  apply  to  leader  propagation. 

The  leader  corona  from  a  350-kV  point  (i.e.,  E  h  —  35  kV/m  x 

10  m)  would  normally  have  (Ref.  22)  a  charge  of  several  uC, 

which  would  thus  draw  a  conduction  current  of  about  10  A  in 
3  x  10  5  mho/m  ionized  air.  This  is  typical  of  measured  leader 
channel  currents  in  long  gaps  (Refs.  22,23,24).  Leader  channel 
velocities  are  typically  observed  to  be  related  to  current  by 
the  empirical  rule  (Ref.  22)  v  (m/s)  ^  2  x  10^  i  (A) .  Thus  a 
10 -A  current  implies  a  velocity  of  2  x  105  m/s,  which  is  con¬ 
sistent  with  the  observed  leader  velocities  on  the  Mike  shot 
(Ref.  15).  Finally,  the  potential  drop  along  a  leader  channel 
for  a  long  spark  is  generally  observed  (Ref.  22)  to  be  of  the 
order  of  (kV/m)  ^  100  i  ,  or  about  30  kV/m  for  a  10-A 
current.  Since  the  EMP  fields  are  calculated  (with  +  50  percent 


36 


uncertainty  in  Jo/oq)  to  ke  °f  this  same  order  of  magnitude 
at  the  earth's  surface  but  will  decrease  with  increasing  ele¬ 
vation,  it  would  not  be  surprising  for  the  Mike  lightning 
strokes  to  stop  in  midair,  which  is  what  may  have  been  observed 
for  four  out  of  the  five  strokes.  It  should  also  be  mentioned 
that  there  are  indications  (Ref.  22)  in  favor  of  even  lower 
potential  gradients  in  very  long  leader  channels.  This  is  sup¬ 
ported  both  by  the  fact  that  cloud- to-ground  lightning  strokes 
commonly  propagate  into  regions  where  Eq  <  20  kV/m  [these  are 
also  convergent  field  geometries  which  are  harder  to  break  down, 
at  least  in  short  gaps,  than  the  divergent  field  geometries 
being  discussed  here  (Refs.  19,22)],  and  by  the  extrapolation 
of  data  on  long  positive-point-to-plane  gaps  (mentioned  earlier) , 
which  implies  that  a  1-m  increase  in  gap  length  at  di  600  m 
only  increases  the  breakdown  potential  by  10  kv. 

Thus,  to  summarize  the  discussion  of  field  strength  criteria 
for  nuclear  lightning  initiation  and  propagation,  it  must  be 
concluded  that  available  data  from  naturally  occurring  ground- 
to-cloud  lightning,  from  lightning  rod  phenomena,  and  from 
breakdown  of  long  gaps  in  air  cannot  be  unambiguously  applied 
to  predict  the  fields  needed  for  creation  of  upward  propagating 
lightning  in  an  EMP  environment.  However,  these  data  clearly 
do  not  rule  out  the  possibility  that  the  observed  Mike  light¬ 
ning  might  be  consistent  with  the  type  of  EMP  theory  presented 
in  this  report.  Moreover,  consideration  of  the  possible  effects 
cf  ionized  air  conductivity  on  the  discharge  process  suggests 
that  lightning  discharges  might  be  very  different  and  perhaps 
also  easier  to  initiate  and  to  propagate  in  an  EMP  source 
region  than  in  un-ionized  air. 

With  regard  to  the  issue  of  lightning  stroke  curvature,  it 
is  not  immediately  clear  whether  the  dominant  electric  field 
guiding  the  direction  of  stroke  propagation  will  be  the  ambient 


EMP  field  or  the  field  due  to  the  stroke  itself.  However, 
since  one  might  expect  any  "self-field"  to  favor  essentially 
vertical  propagation  (whereas  the  strokes  are  actually  curved 
inward,  roughly  concentric  about  the  burst  point),  we  will 
assume  that  the  strokes  are  tending  primarily  to  follow  the 
local  direction  of  the  ambient  EMP  field;  we  may  then  ask 
whether  this  indicated  direction  is  consistent  with  EMP  theory. 

The  predominant  indication  from  the  observed  curvature  is 
that  Eq  >>  E^  throughout  most  of  the  temporal  and  spatial  range 
of  interest.  Apparently  based  upon  Longmire’s  approximate 
analysis  of  Model  2  (ion-dominated  conductivity,  also  see 
Figure  3) ,  it  has  been  suggested  (Ref.  8)  that  this  observa¬ 
tion  implies  a  negligible  electron  conductivity.  Since  at 

the  expected  source  strength  for  a  10-MT  burst  (S  -  3.6  x 

29  30  ® 

10  -  1.1  x  10  ion-pairs/m/s  for  t  »  1-3  ms) ,  electronic 

contributions  would  normally  be  expected  to  dominate  out  to 

r  »  800  m  to  1000  m,  it  has  been  concluded  that  the  lightning 

stroke  curvature  demonstrates  the  need  for  some  anomalous 

mechanisms  for  electron  removal. 

However,  our  present  solutions  for  Model  3  (see  Figure  4), 
with  predominantly  electronic  conductivity  at  closer  ranges 
and  predominantly  ionic  conductivity  farther  out,  clearly 
demonstrate  that  one  can  still  find  regions  where  Eg  >>  Er 
even  when  the  electronic  conductivity  is  not  negligible  through¬ 
out  the  region  of  interest.  In  Figure  7  we  have  shown  the 
angular  dependence  of  Eg  and  E  at  r  =  960  m  and  r  =  1280  m 
as  predicted  by  Model  3.  It  may  be  seen  that  even  as  close 
in  as  r  =  960  m.  Eg  clearly  dominates  over  Er  up  to  eleva¬ 
tions  of  30°  or  more  above  the  ground  plane. 

The  electric  field  structure  reflected  in  Figure  7  shows 
somewhat  less  ir  "ird  curvature  than  is  seen  in  the  photographs 
of  the  Mike  lightning,  primarily  over  the  upper  half  of  the 
more  close-in  observed  strokes.  Thus,  the  ratio  of  ionic  to 
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6  3  kV/m)  . 


electronic  conductivity  during  the  later  stages  of  the  stroke 
propagation  was  probably  somewhat  greater  than  allowed  in 
Model  3,  and  the  ion-dominated  Model  2  may  also  have  become 
applicable  somewhat  earlier  than  we  have  assumed.  However, 
for  the  reasons  given  below,  these  observations  are  not 
inconsistent  with  the  present  theory. 

First,  it  should  be  noted  that  the  ionization  rate  would 
have  fallen  by  about  a  factor  of  3  over  the  t  ~  1-  to  10-ms 
time  frame  during  which  the  lightning  strokes  evolved,  thereby 
increasing  the  ratio  o^/o  .  Furthermore,  it  must  be  realized 
that  the  parameter  values  used  in  the  above  calculations  have 
simply  reflected  our  best  a  priori  estimates.  They  are  also 
subject  to  appreciable  uncertainty.  We  have  already  mentioned 
the  estimated  ±30  to  50  percent  uncertainties  (at  a  given 
water  vapor  content  and  electric  field  strength)  in  the  values 
of  electron  attachment  rate,  electron  and  ion  mobilities,  and 
ion-ion  recombination  rate  constant.  In  addition,  the  value 
of  Oj/o  depends  on  the  dose  rate  (which  is  probably  not  known 
at  any  given  time  to  within  a  factor  of  2)  and  the  water  vapor 
concentration  (which  is  probably  also  uncertain  by  ±30  to 
50  percent). 

When  accounting  for  these  uncertainties  in  the  parameters 

which  influence  the  value  of  oT/c  and  for  the  decline  in  dose 

-L  e 

rate  throughout  the  t  =  1-  to  10-ns  time  frame,  the  details  of 
the  observed  lightning  stroke  curvature  do  not  appear  to  be 
inconsistent  with  the  existing  theory  of  source  region  EMP 
air  chemistry  as  reflected  in  the  present  calculations .  At 
this  time,  no  anomalous  electron-attachment  mechanisms  appear 
to  be  required  to  rationalize  the  observed  curvature. 
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APPENDIX  A.  NUMERICAL  ANALYSIS 


Although  we  have  exact  closed-form  solutions  for  each 
model,  some  numerical  work  is  needed  to  obtain  explicit 
results.  Fortunately,  this  is  not  too  difficult,  and  all 
of  the  results  in  this  report  have  been  obtained  using  a 
Hewlett-Packard  65  programmable  calculator.  The  coefficients 
A^  can  be  calculated  directly  from  Eq.  (10)  ,  and  the  Legendre 
polynomials  and  their  derivatives  can  be  easily  calculated 
from  their  standard  recursion  relations  (Refs.  13,14),  start¬ 
ing  from  Pq(cos  9)  =  1  and  P^(cos  0)  =  cos  0.  The  special 
functions  can  be  calculated  from  their  series  expansions 
(Refs.  10,11) 
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The  function  U(a,b,x)  can  be  expressed  (Ref.  14)  in  terms 
of  M(a,b,x)  and  M(l+a-b,2-b,x) ,  which  may  then  be  evaluated 
from  Eq.  (A3) .  The  gamma  function  may  be  very  accurately  cal¬ 
culated  for  | x |  £  1  from  a  26-term  polynomial  approximation 
(Ref.  11)  for  1/r (x) ,  together  with  the  standard  recursion 


relation  for  |x|  >  1.  Considerable  accuracy  is  sometimes  needed 
for  T (x)  in  order  to  evaluate  the  first  integral  in  Eq.  (18) . 


The  various  integrals  can  be  evaluated  by  substituting 
the  above  series  expansions  for  the  special  functions  and  for¬ 
mally  integrating  term-by-term  to  obtain  a  series  expansion 
for  the  integral.  For  Models  1  and  3  one  has  an  integral  of 
the  form 
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which  may  be  evaluated  by  numerical  summation,  and 
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which  may  also  be  numerically  summed  with  the  exponential 
integral  En(x)  evaluated  by  (Ref.  14) 


Eo(x>  -  le'X 


and 


E1(x) 


.5772156649.  . . 


£n  (x) 


CO 


m=l 


(A6 ) 


44 


Em+l(x)  =  in  (e 
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For  Model  2  one  has  integrals  of  the  form 
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which  may  be  numerically  summed  with  the  incomplete  gamma 
function  Y(a,x)  evaluated  from  (Ref.  14) 


and 
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For  Model  2  one  also  has  integrals  which  are  obtained  by 
expressing  U(a,b,x)  in  terms  of  M(a,b,x)  and  M (1+a-b, 2-b,x) 
and  which  are  of  the  form 

J  dy  M (a ,b,y) yce  2y  = 

x£  (l)^1  JT  (r  (n+a+l)  -  Y  (n+a+1  rx)) 

n=0  n 


which  may  be  evaluated  by  numerical  summation  with  the  gamma 
and  incomplete  gamma  functions  evaluated  as  described  above. 
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For  each  of  the  three  models,  ^(r)  and  d<fJ,(r)/dr  are 
calculated  separately  from  their  respective  closed-form 
expressions.  Thus,  a  useful  cross-check  of  the  numerical 
results  is  obtained  by  comparing  the  calculated  values  of 
dtji^rj/dr  with  values  obtained  by  finite  differences  from  the 
calculated  values  of  ^(r)  at  nearby  range  points. 

The  numerical  results  for  4>^(r)  and  d(j>^(r)/dr  are  dis¬ 
played  in  dimensionless  form  in  Figures  A1  through  A6.  It  may 
be  seen  that  in  some  cases  the  expansion  in  Legendre  polynomials 
will  not  converge  particularly  rapidly,  so  that  the  accuracy 
with  which  the  electric  fields  are  calculated  may  sometimes  be 
only  ~10  percent.  In  an  attempt  to  improve  convergence  with  a 
minimum  of  computational  effort,  values  for  4>^(r)  have  been 
estimated  for  the  highest  SL  by  extrapolation  of  the  values  for 
lower  l.  In  each  figure,  the  absolute  value  is  plotted,  with 
the  sign  denoted  by  the  (+)  or  (-)  symbol. 
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igure  A2.  Dimensionless  results  for  d/dr  <j>.  (r) 
for  Model  1. 
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Figure  A5.  Dimensionless  results  for  <f>{  (r)  for 
Model  3. 
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